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The realization of a scalable quantum information processor has emerged over the past decade 
as one of the central challenges at the interface of fundamental science and engineering. Much 
progress has been made towards this goal. Indeed, quantum operations have been demonstrated 
on several trapped ion qubits, and other solid-state systems are approaching similar levels of con- 
trol. Extending these techniques to achieve fault-tolerant operations in larger systems with more 
qubits remains an extremely challenging goal, in part, due to the substantial technical complexity 
of current implementations. Here, we propose and analyze an architecture for a scalable, solid-state 
quantum information processor capable of operating at or near room temperature. The architecture 
is applicable to realistic conditions, which include disorder and relevant decoherence mechanisms, 
and includes a hierarchy of control at successive length scales. Our approach is based upon re- 
cent experimental advances involving Nitrogen- Vacancy color centers in diamond and will provide 
fundamental insights into the physics of non-equilibrium many-body quantum systems. Addition- 
ally, the proposed architecture may greatly alleviate the stringent constraints, currently limiting the 
realization of scalable quantum processors. 



Nitrogen- Vacancy (NV) color centers in diamond stand 
out among other promising qubit implementations ^-[7] 
in that their electronic spins can be individually polar- 
ized, manipulated and optically detected under room- 
temperature conditions. Recent advances involving the 
quantum manipulation of such crystal defects have al- 
lowed researchers to achieve sub-diffraction limited res- 
olution and dipole-coupling mediated entanglement be- 
tween neighboring NV electronic spins [BHTU]. Despite 
such substantial developments in this and other experi- 
mental systems, it remains unclear whether these pieces 
can be combined into a scalable quantum information 
processor (QIP) capable of operating under ambient, 
room temperature conditions. 

In what follows, we describe and analyze a feasible 
architecture for a diamond-based quantum information 
processor. Our approach makes use of an array of single 
NV centers, created through implantation of ions and 
subsequent annealing [TTl [20]. Each NV center consti- 
tutes an individual quantum register containing a nu- 
clear spin and a localized electronic spin. The nuclear 
spin, which has a long coherence time, serves as the 
memory qubit, storing quantum information, while the 
electronic spin will be used to initialize, read out, and 
mediate coupling between nuclear spins of adjacent reg- 
isters. Magnetic dipole interactions allow for coherent 
coupling between NV centers spatially separated by tens 
of nanometers. While in principle, a perfect array of NV 
centers would enable scalable quantum information pro- 
cessing, in practice, the finite creation efficiency of such 
centers, along with the requirements for parallelism, ne- 
cessitate the coupling of registers separated by signifi- 
cantly larger distances. To overcome this challenge, we 



show that the coupling between NV centers can be me- 
diated by an optically un-addressable "dark" spin chain 
data bus (DSCB). For concreteness, within our architec- 
ture, we will consider the specific implementation of such 
a DSCB by utilizing implanted Nitrogen impurities (P3 
centers) with spin 1/2, as shown in Fig. la. [51 121j. 

The NV Qubit Register 

Single NV registers contain a spin triplet electronic 
ground state [S = 1) and can be optically pumped and 
initialized to the |0)g spin state, which has no magnetic 
dipole coupling with other NV registers or impurities. Af- 
ter optical initialization, the electronic spin of each regis- 
ter remains in the |0)g state, unless coherently transferred 
to the |l)g state by a resonant microwave (MW) pulse, 
as depicted in Fig. la [T0Hl3l . The nuclear spin associ- 
ated with Nitrogen atoms (/ = 1/2 for ^^N) possesses 
an extremely long coherence time and will serve as the 
memory qubit in our system |23| : manipulation of the 
nuclear spin is accomplished with RF pulses [22]. The 
Hamiltonian governing the electronic and nuclear spin of 
the NV register is 

= Ao^^ ^ ^^^^^ ^ ^^^Qj^ ^ Agj^_^ 

with zero-field splitting Aq — 2.87GHz, electronic spin 
gyromagnetic ratio /ie = — 2.8MHz/Gauss, nuclear spin 
gyromagnetic ratio /i„ = —0.43 kHz/Gauss, and hyper- 
fine coupling A = 3.0 MHz (TU]. The application of a 
magnetic field along the NV-axis {z) ensures full address- 
ability of the two-qubit system, resulting in the energy 
levels shown in Fig. la. A universal set of two qubit quan- 
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FIG. 1: Schematic representation of individual NV registers 
within bulk diamond, (a) Each NV register contains a nu- 
clear spin 1 = 1/2 (yellow), providing quantum memory, and 
an electronic spin 5=1 (green). Dark spins (black) repre- 
sent elements of an optically un-addressable spin chain which 
coherently couples spatially separated NV registers. The NV 
level structure (in a high B field) is shown. A resonant mi- 
crowave {^Mw) pulse coherently transfers the electronic spin 
of the register from \Q)^ to 11)^; subsequent manipulation of 
the nuclear spin is accomplished through an RF pulse (f2_RF). 
The far detuned | — 1)^ state can be neglected to create an 
effective two-qubit register. However, the full three level NV 
structure will be utilized in horizontal DSCB mediated co- 
herent coupling of NV registers, (b) A universal set of two- 
qubit gates can easily be achieved with only MW and RF 
controls [22j. Electronic spin manipulation can be accom- 
plished with a MW field, where t represents the duration of 
the MW pulse. By exploiting the hyperfine coupling between 
the electronic and nuclear spin, one can achieve controlled- 
NOT operations conditioned on either spin. In particular, a 
CeNOTn gate can be accomplished by utilizing a RF vr-pulse, 
which fiips the nuclear spin conditioned on the electronic spin 
being in |l)e. Similarly, a CnNOTe gate can be accomplished 
by utilizing the hyperfine interaction to generate a controlled- 
phase (CP) gate, where r represents the duration of the wait 
time required to achieve such a hyperfine-driven CP gate. 
Performed between two single-qubit Hadamard gates (7r/2- 
pulses) on the electronic spin, such a CP gate generates the 
desired C„NOTe gate. Finally, combining the CeNOT„ and 
CnNOTe gates allows for the execution of a SWAP gate. 



turn operations can easily be achieved with only MW and 
RF controls, as shown in Fig. lb [22) . 

Furthermore, it is possible to selectively readout the 
state of the NV register; for example, to readout the nu- 
clear qubit of a register, we apply a C„NOTe gate to 
couple the electronic and nuclear spins, thereby allowing 
for readout of the electronic spin based on fluorescence 
detection. In the case where NV registers are separated 
by sub-optical- wavelength distances, the readout of regis- 
ters will be complicated by the strong fluorescence back- 
ground from neighboring NV centers. To suppress this 
background fluorescence, a red donut beam can be used, 
with its minimum located at the particular NV center 
being read out [15]. Thus, while the fluorescence signal 
from the NV register located at the minimum persists, 
the remaining illuminated registers will be dominated by 
the stimulated emission induced by the red donut beam. 
In addition to suppressing the background noise, the red 
donut beam can also suppress the nuclear decoherence of 
the remaining NV registers, by reducing the amount of 



time these registers spend in the excited electronic state. 
After each round of fluorescence detection, the electronic 
spin is polarized to the |0)g state, while the 1^ component 
of the nuclear spin, a quantum non-demolition observ- 
able, remains unchanged [51|. Therefore, it is possible 
to repeat this readout procedure multiple times in or- 
der to improve the readout fidelity [T71 [TB]. A strong 
magnetic field ~ 1 Tesla along the NV axis should 
be used to decouple the electronic and nuclear spins in 
order to achieve high fidelity single shot readout of NV 
registers [H]. In addition to sub- wavelength readout, 
optical donut beams also introduce the possibility of se- 
lectively manipulating individual NV registers with sub- 
wavelength resolution. While un-illuminated NV centers 
may respond to a resonant MW pulse, illuminated reg- 
isters undergo a strong optical cycling transition which 
suppresses their response to microwave pulses due to the 
quantum Zeno effect [25l - f27] . 



Approach to Scalable Architecture 

One of the key requirements for fault-tolerant quan- 
tum computation is the ability to perform parallel gate 
operations. In our approach, this is achieved by con- 
sidering a hierarchy of controllability. The lowest level 
of the hierarchy consists of an individual optically ad- 
dressable plaquette with horizontal and vertical spatial 
dimensions ^ 100 — 500nm, containing a single compu- 
tational NV register, as shown in Fig. 2a. The plaque- 
tte dimensions are chosen such that register control and 
readout can be achieved using conventional far-field or 
sub- wavelength optical techniques [TO l (M l ITS l l23 t [27 l [28] . 
The second level, termed a super-plaquette (~ 10/xm 
X lO^m) , consists of a lattice of plaquettes whose com- 
putational registers are coupled through DSCBs. At the 
highest level of the hierarchy, we consider an array of 
super-plaquettes, where individual super-plaquettes are 
controlled by confined microwave fields 29J. In partic- 
ular, micro-solenoids can confine fields to within super- 
plaquettes, allowing for parallel operations at the super- 
plaquette level. For example, as depicted in Fig. 2, in- 
dependent microwave pulses can allow for simultaneous 
operations on the electronic spins of all computational 
NV registers within all super-plaquettes. In order to 
control registers at the super-plaquette boundaries, we 
define a dual super-plaquette lattice (Fig. 2a). Localized 
microwave fields within such a dual lattice can provide a 
smooth transition between the boundaries of neighboring 
super-plaquettes. 

Taking advantage of the separation of length scales 
inherent to optical control and microwave confinement 
provides a mechanism to achieve parallelism; indeed, the 
hierarchical control of plaquettes, super-plaquettes, and 
super-plaquette arrays allows for simultaneous single- 
and two-qubit gate operations, which are fundamental 
to fault-tolerant computation. One of the key difference 
in the currently proposed architecture as compared to 
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FIG. 2: The architecture for a room-temperature solid-state 
quantum computer, (a) A two-dimensional hierarchical lattice 
allowing for length-scale based control, which enables fully 
parallel operations. At the lowest level, individual plaque- 
ttes are outlined in grey and each contains a single computa- 
tional NV register. At the second level of hierarchy, a super- 
plaquette, outlined in white, encompasses a lattice of pla- 
quettes; each super-plaquette is separately manipulated by 
micro-solenoid confined microwave fields. In order to allow 
for quantum information transfer across boundaries of super- 
plaquettes, there exists a dual super-plaquette lattice outlined 
in red. (b) The schematic NV register implantation within a 
super-plaquette. Two rows of individual plaquettes within a 
super-plaquette are shown. NV registers, consisting of an elec- 
tronic (green) and nuclear (yellow) spin are depicted within a 
staggered up-sloping array which is row-repetitive. Individual 
rows within a single plaquette are specified by an integer n 
with n — 1 being the bottom row and n — A'l being the top 
row. To achieve a staggered structure, we specify a unique 
implantation row within each plaquette wherein single impu- 
rities are implanted and subsequently annealed. For a given 
row of plaquettes, the implantation row corresponding to the 
left-most plaquette is n = 1, while the plaquette to the imme- 
diate right has implantation row n = 2; this pattern continues 
until the final plaquette in a given row, which by construc- 
tion, has the highest implantation row number. The implan- 
tation process is repeated for each row of plaquettes within the 
super-plaquette and creates an array of NV registers, which 
each occupy a unique row in the super-plaquette. Since each 
NV register occupies a unique row within the super-plaquette, 
the magnetic field gradient in the y direction allows for indi- 
vidual spectroscopic addressing of single registers. Coherent 
coupling of spatially separated NV registers in adjacent pla- 
quettes is mediated by a dark spin chain data bus (DSCB) 
and is schematically represented by the curved line connect- 
ing individual registers (SI Text). The second implantation 
step corresponds to the creation of these horizontal and ver- 
tical dark spin chains. 



previous proposals [U [3U] is that the design here does 
not rely on optically resolved transitions, which are only 
accessible at cryogenic temperatures. 

The required 2D array of NV centers can be created 
via a two-step implantation process and the selective ma- 
nipulation of individual registers within such an array is 
enabled by the application of a spatially dependent ex- 
ternal magnetic field {y) = + B^^. The ID 
magnetic field gradient is sufficiently strong to allow for 
spectroscopic microwave addressing of individual NV reg- 
isters, each of which occupies a unique row in the super- 
plaquette, as shown in Fig. 2b. 



Results and Discussion 

Dark Spin Chain Data Bus 

To coherently couple two spatially separated NV cen- 
ters, we consider two distinct approaches. First, we con- 
sider an approach, which is appropriate for spin-state 
transfer along the direction of the magnetic field gradi- 
ent, in which individual addressing of spins is possible. 
This allows for an adiabatic sequential SWAP between 
neighboring qubits and, consequently, between the ends 
of the chain. Alternatively, in the situation where indi- 
vidual addressing of spins is not possible (i.e. direction 
transverse to the field gradient), we show that global con- 
trol pulses achieve effective Hamiltonian evolution, which 
enables quantum state transfer through the spin chain. 
In both cases, we show that perfect state transfer and re- 
mote coupling gates are possible even when the interme- 
diate spin chain is completely unpolarized (infinite spin 
temperature). 

We begin by analyzing the adiabatic sequential SWAP 
in a spin-1/2 chain. This approach is suitable to couple 
registers in plaquettes that are vertically adjacent, rely- 
ing upon the individual addressability of qubits and uti- 
lizing the magnetic dipole coupling between spin-chain 
elements. As shown in the SI Text, under the secular 
approximation, the magnetic dipole coupling between a 
pair of neighboring spins can be reduced to Ising form 



Hint — ^nSlSl + 



(^0 + 5i)Sl 



(2) 



where k is the relevant component of the dipole tensor. 
Wo captures the electronic Zeeman energy, and 5i charac- 
terizes both the hyperfine term (nuclear spin dependent) 
and the magnetic field gradient (SI Text). From the Ising 
Hamiltonian, an XX interaction between qubits can be 
distilled by driving with Hdrive — J2i=i i 2f^i5'^ cos[(a;o-l- 
(5,;)t], leading to (under the rotating wave approxima- 
tion, in the rotating frame, and in a rotated basis with 
(x,?/,z) {z,-y,x)) 

H,nt = n{s+s^ + s^s+) + n,sl + n^sl (a) 

The spin-flip process in Hint is highly suppressed in the 
limit of |rii— r22| ^ K, while the same process is dominant 
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in the case of jfii — ^ Hence, by slowly ramp- 
ing the Rabi frequencies Vli and through one another, 
adiabatic SWAP of the quantum states of the two impu- 
rities can be achieved through rapid adiabatic passage, as 
shown in Fig. 3a. Generalizing to arbitrary length spin 
chains yields H^* = ^(St S- S+^^)+j:^ n,Sl, 
whereby the sequential adiabatic SWAP of quantum 
states along the spin chain can be achieved by succes- 
sively tuning individual Rabi frequencies across one an- 
other. During the adiabatic SWAP of a single pair of 
spins, higher order interactions, such as those resulting 
from next-to-nearest neighbors, will be suppressed due 
to the differences in Rabi frequencies. By including the 
magnetic dipole coupling between the electronic spin of 
the NV register and the spin chain quantum channel, we 
arrive at an effective mixed spin chain with the DSCB 
connecting the two electronic spins of the vertically sep- 
arated NV registers. The specific procedure resulting in 
adiabatic sequential SWAP mediated coupling between 
NV registers is depicted in Fig. 3c. 

Crucially, such an adiabatic sequential SWAP is robust 
against variations in the coupling strength k, which can 
be induced by the imprecise implantation of impurities 
that form the spin-1/2 chain; in particular, even for the 
case of varying K^.i+i, perfect adiabatic SWAP occurs so 
long as the rate at which fli and ili+i are ramped through 
one another is sufficiently small. Within the proposed 
architecture, the impurities forming the horizontal spin 
chain will not induce operational errors during the verti- 
cal adiabatic sequential SWAP since the design principle 
allows for selective spin echoing (SI Text). 

Next, we consider a second method, termed free 
fermion state transfer (FFST) developed in [31], to coher- 
ently couple NV registers. In contrast to the adiabatic 
sequential SWAP, the method utilizes only global con- 
trol over impurities and effective Hamiltonian evolution. 
The relaxation of the requirement of individual control 
over elements of the dark spin chain renders this second 
method applicable for coherent coupling between NV reg- 
isters in horizontally adjacent plaquettes, transverse to 
the direction of the field gradient. In particular, the pro- 
tocol achieves coherent coupling through an unpolarized, 
infinite temperature spin chain, employing purely Hamil- 
tonian evolution under 



Magnetic Field Gradient 



HfFST = g{S^Vi^l + ^NV2^N + h.C. 



N-1 



(4) 



as shown in Fig. 3b. This Hamiltonian, obtained in anal- 
ogy to Eq. 3, results in coherent interactions between 
NV centers, which is best understood via an analogy to 
eigenmode tunneling in a many-body system. Specifi- 
cally, the spin chain described by Hffst can viewed as a 
system of non- interacting fermions. As described in |31j . 
by tuning the NV centers into resonance with a single 
fermionic eigenmode, an effective three-state system can 
be realized. Mediated by this fermionic eigenmode, the 
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FIG. 3: Dark spin chain data bus (DSCB) mediated coher- 
ent coupling of spatially separated NV registers, which does 
not require spin chain initiahzation. (a) Adiabatic sequential 
SWAP along the vertical direction, parallel to the magnetic 
field gradient. Individual addressing of impurities, enabled 
by the field gradient, allows for a slow ramping of the Rabi 
frequencies Hi and Qj through one another; this achieves adi- 
abatic SWAP of the quantum states of the two impurities 
through rapid adiabatic passage. Thus, sequential adiabatic 
SWAP of quantum states along the spin chain can be achieved 
by successively tuning individual Rabi frequencies across one 
another, (b) Free fermion state transfer in the horizontal di- 
rection, transverse to the magnetic field gradient. The cou- 
pling strength between the end qubits and the spin chain is g, 
while the interchain coupling strength is k. Schematic repre- 
sentation of the level structure of the NV electronic spin and 
a dark impurity spin. Controlling the NV-impurity coupling g 
is an essential component of FFST and occurs by driving the 
NV in two-photon resonance, with Rabi frequency Q and de- 
tuning A. (c) Schematic circuit diagram outlining the proto- 
col to achieve coherent coupling between the nuclear memory 
qubits of spatially separated NV registers. First, the nuclear 
and electronic qubits of a single register are swapped. Next, 
the electronic qubits of the two NV centers to be coupled are 
swapped via the DSCB. Finally, a two-qubit gate between the 
electronic and nuclear spin of the second register is performed 
before the memory qubit is returned to the nuclear spin of the 
original NV center. 
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electronic states of two remote NV centers can be co- 
herently swapped, leading to an analogous protocol for 
remote register coupling as shown in Fig. 3c. Crucially, 
such a SWAP gate is insensitive to the polarization of the 
intermediate dark spins and high-fidelity quantum state 
transfer can be achieved, provided that the fermionic 
mode is delocalized and that the coupling, g, of the NV 
qubit to the spin chain is controllable. As detailed in the 
Materials and Methods, by utilizing the three-level NV 
ground-state structure (Fig. 3b), it is possible to fully 
control the NV-chain coupling. This tunability also en- 
sures that FFST is fundamentally robust to experimen- 
tally relevant coupling-strength disorder, which could be 
induced by implantation imprecision. Indeed, by sepa- 
rately tuning the NV-chain coupling on either side of the 
DSCB, it is possible to compensate for both disorder- 
induced asymmetry in the fermionic eigenmode as well as 
altered statistics of the eigenenergies [3Tti33] . The partic- 
ular implementation of FFST within the NV architecture 
is further described in the SI Text. 



Implementation, Operational Errors and Gate Fidelities 

The specific implementation of the DSCB can be 
achieved with implanted Nitrogen impurity ions. Dipole 
coupling between neighboring Nitrogen electronic spins 
forms the DSCB, while dipole coupling between the NV 
and Nitrogen electronic spins forms the qubit-DSCB in- 
teraction; non-secular terms of this magnetic dipole cou- 
pling are highly suppressed due to the spatially depen- 
dent external magnetic field Bz{y), resulting in the ef- 
fective interaction found in Eq. 2. In addition, the Ni- 
trogen impurities possess a strong hyperfine coupling, 
the principal axis of which can take on four possible 
orientations due to tetrahedral symmetry [341 - I36] . Dy- 
namic Jahn- Teller (JT) reorientation of the Nitrogen im- 
purity's hyperfine principal axis results in two particu- 
lar considerations: 1) the addressing of additional JT 
frequencies yielding a denser super-plaquette frequency 
spectrum and 2) the JT-governed spin-lattice relaxation 
(SLR) time (SI Text). Since is characterized by 
an Arrhenius rate equation |35| at ambient temperatures, 
a combination of a static electric field and slight cooling 
by w 50K allows for a substantial extension of the relax- 
ation time to ~ Is; hence, in the following consideration 
of operational errors, we will assume that we are lim- 
ited by T(^^ , the spin-lattice relaxation time of the NV 
center. 

We now consider various imperfections, which may in- 
troduce operational errors. In particular, we consider 
the errors associated with: 1) sequential SWAP medi- 
ated coupling between vertically adjacent registers and 
2) FFST between horizontally adjacent registers. We be- 
gin by discussing the analytic error estimate associated 
with each method, after which, we summarize the results 
of full numerical simulations. 

First, we consider the accumulated infidelity associated 
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FIG. 4: Numerical simulation of the DSCB fidelity, (a) The 
operational infideUty associated with the adiabatic sequen- 
tial SWAP for A'' = 18. The simulations account for the 
Jahn- Teller orientation of Nitrogen impurities and utilize the 
optimized adiabatic ramp profile [37] ■ Simulations utilize 
an optimized coupling strength of 8.71kHz (IS.lnm spacing). 
Full numerical integration of the time dependent Schrodinger 
equation produces infidelity contour plots as a function of to- 
tal SWAP time and T^^. (b) Numerical simulations of the 
operational infidelity associated with FFST for N = 7. Non- 
nearest neighbor interactions are assumed to be refocused 
through dynamic decoupling as described in the SI Text. 
Simulations, which utilize an optimized coupling strength of 
12.6kHz (16nm spacing), are based upon a full diagonalization 
and also account for the Jahn- Teller orientation of Nitrogen 
impurities. Infidelity contour plots are again shown as a func- 
tion of total SWAP time and Tf^ . 



with the adiabatic sequential SWAP (SI Text), 



P. 



ss 



N (p^ff + pad^a + PMp + PtI + PT2 ) • (5) 



The first term, p^s 



off '~ (s^) ' represents off-resonant 
excitations induced by microwave manipulations with 
Rabi frequency 51^. Here, Ag characterizes the gradient- 
induced splitting achieved within the super-plaquette fre- 
quency spectrum. The second term, Padia, corresponds 
to the non-adiabatic correction resulting from an opti- 
mized adiabatic ramp profile [371 - 139] . The third term, 

Pdip ^ (s?^) ' directly obtained from Eq. 3 and cor- 
responds to additional off-resonant errors. The fourth 
error term, p^f corresponds to the depolarization error 
induced by the finite NV Ti time, while the final error 
term, p^f corresponds to the infidelity induced by de- 
phasing. Since each error term is considered within the 
context of a single adiabatic SWAP, the total error con- 
tains an additional factor of N, representing the chain 
length, which is plaquette size dependent (e.g. TV « 5 for 
lOOnm and TV w 20 for 500nm). 

We can similarly consider the accumulated infidelity 
associated with FFST (SI Text), 



„FFST 



„FFST 
Poff 



Pf^ 



- Pg + Pti 



„FFST 
-PT2 



(6) 



In direct analogy to pf^- 1 the first term in pf^^^ corre- 
sponds to the excitation of an NV register by off-resonant 
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microwave fields. The second term, pfermi, corresponds 
to the undesired couphng with off-resonant fermionic 
modes. Since the coupling strength is characterized by 
g/VN ^T], while the splitting of the eigenenergy spec- 
trum ~ k/N, such an off-resonant error induces an infi- 
delity ^ i ^i^f^ )^ ■ The third error term, pg, results from 
the protocol designed to control, g, the NV-chain cou- 
pling (see Materials and Methods for details). Finally, 
directly analogous to pfrn the fourth and fifth terms 
correspond to errors induced by the operational time, 
tpFST, which causes both depolarization and dephasing. 

Finally, we perform numerical simulations, taking into 
account the Nitrogen JT frequencies, to characterize the 
infidelity of both the adiabatic sequential SWAP and 
FFST within the NV architecture, as depicted in Fig. 4. 
The results of these calculations are in excellent agree- 
ment with the above theoretical predictions (SI Text). 
In particular, these simulations reveal that, for suffi- 
ciently long ~ 100ms, operational infidelities in 
both DSCB methods can be kept below 10~^. 

These simulations clearly show that the Ti time of the 
NV electronic spin is of critical importance in obtaining 
high-fidelity quantum operations. While at room tem- 
perature Ti appears to vary depending on the particular 
sample and on the specific properties of the local NV en- 
vironment, such as strain, values on the order of 10ms 
are generally obtained [231 IMj ■ However, the spin-lattice 
relaxation mechanism governing Ti is most likely related 
to an Orbach process [101 EI] , which is strongly temper- 
ature dependent. In such a case, modest cooling of the 
sample by « 50K, is likely to extend TI by more than an 
order of magnitude, thereby making high fidelity gates 
possible. 

Given that such numerical estimations suggest the pos- 
sibility of achieving high fidelity two-qubit operations be- 
tween remote NV registers, the proposed architecture 
seems well suited to the implementation of topological 
quantum error correction [42H45| . Recent progress in op- 
timizing the 2D nearest-neighbor surface code has yielded 
an error threshold of e « 1.4% [IS], which is above the 
estimated infidelity corresponding to both the adiabatic 
sequential SWAP and FFST; thus, in principle, imple- 
mentation of the 2D surface code can allow for successful 
topological quantum error correction, and hence, fault 
tolerant quantum computation [47] . 

Conclusions and Outlook 

The above considerations indicate the feasibility of ex- 
perimentally realizing a solid-state quantum computer 
capable of operating under ambient conditions at or near 
room temperature. We emphasize that a majority of the 
elements required for the realization of individual qubits 
in our architecture have already been recently demon- 
strated. In our approach, these techniques are supple- 
mented by both a new mechanism for remote register 
coupling between NV centers as well as a hierarchical de- 



sign principle, which facilitates scalability. The remote 
coupling mechanisms discussed can naturally be imple- 
mented via Nitrogen ion implantation in ultra-pure di- 
amond crystals and are robust to realistic imperfections 
and disorder [ST] . 

While the implementation and integration of the vari- 
ous proposed elements still require substantial advances 
in areas ranging from quantum control to materials sci- 
ence, a feasible approach to room temperature quantum 
information processing can greatly alleviate the strin- 
gent requirements associated with cryogenic tempera- 
tures, thereby making the realization of a scalable quan- 
tum computer significantly more practical. 

The present work opens a number of new directions 
which can subsequently be explored. In particular, while 
we have considered the direct errors associated with 
DSCB mediated coupling, it is instructive to note that 
the fidelity of such quantum gates can often be signif- 
icantly improved using techniques from optimal control 
theory [iSUSOj . For example, such methods of optimal 
control, while negating the detrimental effects of deco- 
herence, can also simultaneously allow for the implemen- 
tation of high-fidelity gates despite both frequency and 
coupling disorder as induced by ion implantation errors. 
Indeed, the ability to precisely guide the quantum evolu- 
tion via optimal control, even when the system complex- 
ity is exacerbated by environmental coupling, provides 
an alternative solution to improve single and two-qubit 
gate fidelities fST. In addition, it is well known that the 
local strain field surrounding each NV center can signif- 
icantly alter the register's properties; hence, through a 
detailed understanding of electric field induced strain, it 
may be possible to improve the coherence properties of 
the qubit. Beyond these specific applications, a number 
of scientific avenues can be explored, including for exam- 
ple, understanding and controlling the non-equilibrium 
dynamics of disordered spin systems. 

Materials and Methods 

Controlling Qubit-Chain Coupling in the NV 
Architecture — To achieve an effective Hamiltonian of 
the form given by Eq. 4, it is essential to control the 
coupling strength between the NV register and the neigh- 
boring impurity. Here, we utilize the three levels of the 
NV electronic spin [55] to effectively control g, as shown 
in Fig. 3b, whereby the Hamiltonian (under microwave 
driving) can be written as 

H= -A(|l)(l| + |-l)(-l|) 

- r!(|0)(l| + |0)(-1| + h.c.) - nNS^ + 4«:5f , 

(7) 

where represents the Rabi frequency on the NV reg- 
ister, A represents the associated detuning, and 
represents the Rabi frequency on the Nitrogen impu- 
rity. In this case, since the NV two-photon detuning 
is zero, it is convenient to define bright and dark states. 
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\B) 



72 



and |D) "'"^ ; further, in the result- 



ing two-level picture, the associated dressed states are 
1+) « \B) + 4^|0) and I-) « |0) - ^\B), in the limit 
SI <C A. Hence, rewriting the Hamiltonian in this limit 
yields 



(A+^)|+)(+| 



H= -^D){D\ 

+ 2k{\B){D\ + \D){B\){\+)m{ 



1 



-)n{+\), 



(8) 

where |±)jv = '"'"^"^•'-^'^ correspond to the two S^- 
eigenstates of the Nitrogen impurity. The coupling term 
can be further re-expressed as 



2k{{\+){D\ + \D){+\) - 



i\-){D\ + \D){-\ 



'N 



-N 



-N 



){+n\). 



(9) 

Thus, by working within the NV subspacc {\D), |— )}, it 
is possible to completely control the coupling between the 
NV register and Nitrogen impurity, g ~ k£, by tuning 
the Rabi frequency and detuning. It is possible to work in 
the required two-state subspace by ensuring that k <C A 
and hence, that the |-|-) state remains unpopulated, with 
corresponding off-resonant error k^/A^. 

Furthermore, we evince a possible scheme to coher- 



ently map the quantum information that is stored in 
the nuclear memory into the desired electronic subspace. 
For example, consider mapping |0) (a| f) + P\ 4-)) to 
(a|— ) + (3\D)) (8) I t)) where the first (tensor) factor corre- 
sponds to the electronic state and the second corresponds 
to the nuclear state of a single NV. The proposed map- 
ping can be achieved in a two-step process. First, by 
simultaneously performing a tt— pulse on the transitions 

|o) ® I ;) ^ I 1) ® I ;) and |o) ® I ;) ^ |i) ® I ;) with 

oppositely signed Rabi frequencies, one can map |0)(8) | 4-) 
to \D) I J,). Next, one utilizes an RF pulse to flip the 
nuclear spin, which yields \D) (E) \ i) ^ \D) t)- Fi- 
nally, turning Q on in an adiabatic fashion ensures that 
the state preparation populates only \D) and |— ), thereby 
mapping the quantum information into the desired elec- 
tronic subspace. 
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Supporting Information (SI Text) 

Realization of Dark Spin Chain Data Bus with 
Nitrogen Impurities 

Nitrogen Impurity Hamiltonian 

We consider a dipole-dipole coupled chain of nitrogen im- 
purities which forms the basis for the DSCB that couples 
remote NV registers. The Hamiltonian for a single Ni- 
trogen impurity is given by 



Hn = gefJ-sB ■ S - gNtJ-NB ■ I 



(SI) 



where Ay = -159.7MHz and = -113.8MHz are the 
hyperfine constants corresponding to the primed axes, 
which is chosen with the Jahn- Teller ( JT) axis as z' jSl| . 
As in recent reports |S2j . we consider the application of 
a magnetic field B of strength B^ {'^ IT), as defined in 
the main text, throughout the 2D array along the NV 
axis (consider e.g. NV \\ (1,1,1)), which defines the 
quantization axis z. To model the additional Jahn- Teller 
frequencies applied in the super-plaquettes, we calculate 
the impurity chain Hamiltonian between two nitrogens 
with differing Jahn- Teller axes (assuming a (1,1,1) cut 
diamond, as shown in Fig. SI), and with the JT axis of 
A'^i parallel to the NV axis (note that the non-parallel 
JT axes are equivalent to one another relative to the 
NV axis). We explicitly represent the hyperfine inter- 
action within a standard Cartesian basis by choosing the 
Jahn- Teller axis of the first nitrogen iVi as e^/ || (1,1,1) 
and subsequently ex' \\ (2,-1,-1) and Cyi \\ (0,1,-1), 
and of the second nitrogen N2 as Cz' || (1,1,-1) and 
subsequently e^' \\ (-2,1,-1) and Cyi \\ (0,1,1) (the 
generalization to other pairs of JT axes follows after). 
Here, the hyperfine interaction of TVi then takes the form 
Hhf = A\\SzIz + Ai_{SxIx + Syly), wfillc the hyperfine 
term of the TVs takes the form i^y (^ -I- ^ ^)(i -|- 

2i_i_?ii^_i_ A , n( -2S, I 2S„ 4S^w -2J, , 2Iy 4/^ ^ , 

'2S, 2S,_-)(-2i^ _ important to note that the 



V3 

electronic Zeeman term in the impurity Hamiltonian pro- 
vides a shift ~ lOGHz and that the nuclear Zeeman term 
provides a shift ~ IMHz. Thus, to good approximation, 
we can disregard electronic spin flip terms proportional 
to Sx and Sy since the hyperfine constants are two orders 
of magnitude smaller than the electronic Zeeman energy, 
yielding Hhf ~ AnSzIz for A^i and 



H 



HF 



Szlzi 



.Aw 



8A| 



) + Sz{ 



All 



A_^^^2h 



2L. 



3^3 3V3 V6 V2 
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(S2) 

for N2, where {a, 13,-/) = (-7.2, -12.5, -118.9)MHz. For 
iVs, the hyperfine subspace is separated into two mani- 
folds of nuclear spin 1/2 and —1/2; to see that the nuclear 



spin-flip terms are highly off-resonant, it is easiest to work 
in a rotated basis with z = , _ ^ _ =^{ax+(3y+jz), x = 



= (—72; + az), and y 



li.l+\)y + z), wherein Hhf ~ ^^a^ + (3^ + j^Szh- In 
the rotated basis, nuclear spin flips are captured by the 



Zeeman term —gjyiiNBoI^ 

a J J 



-.g7VMAr-Bo( 



; -ij. I , =Iv)- Thus, for a given 

electronic spin, the nuclear spin subspaces are separated 
by ^/o^ + (P~^^ ~ lOOMHz, while the strength of 
the nuclear spin flip terms proportional to and ly 
^ lOOkHz and can hence be neglected. These secular ap- 
proximations yield a nearest neighbor Hamiltonian given 
by (returning to the original basis for simplicity of nota- 
tion) 



Hmi = gef^BBo{Sl + si) — gNl^NBo{ll + II) 



+ A\^S'j'^+Sili{ 



A« 8A_L 



(S3) 



Mo 



+ i^,9lAiSlSl) 



where we have additionally neglected electronic spin 
flip terms originating from dipole-dipole interactions be- 
tween the impurities since their hyperfine terms vary 
by ~ lOMHz while the dipole coupling strength is only 
~ lOkHz for 20nm separation (neglecting the field gradi- 
ent). Thus, in the case where neighboring impurities have 
differing Jahn- Teller axes, the electronic spin flip terms 
can always be neglected to give an Ising interaction. In 
the case of parallel Jahn- Teller axes, when the neigh- 
boring nitrogen nuclear spins are different, the combined 
hyperfine term takes the form ±{Sl-Sl){^ + ^). The 
dipole-dipole terms which correspond to electronic spin 
hipping will attempt to couple | fi , 4-2) and | ii , t2) ; how- 
ever, this coupling is again highly suppressed since these 
states are separated by (^ -f ^4^) ^ lOOMHz, while the 
dipole coupling strength is ~ lOkHz. Furthermore, the 
states I tijt2) and | ii,i2) are separated from the others 
by the electronic Zeeman energy, ensuring that dipole- 
dipole induced spin flips will again be highly off- resonant. 
Finally, when the Jahn- Teller axes are parallel and the 
nuclear spins are also identical, these dipolar spin-flip 
terms are suppressed by the external magnetic field gra- 
dient. It is important to note that for each Nitrogen 
spin, since the hyperfine term can take on four possible 
values, all four frequencies, corresponding to ± \-^\\ 

and ± ^ ( 4f + ) are applied in order to address the 
impurity (wq as defined in Eq. 2 of the main text). 



Dynamic Jahn-Teller Re- orientation 

As mentioned in the main text, it is essential to con- 
sider the intrinsic properties of the dark Nitrogen im- 
purity spin chain; in particular, the room temperature 
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NV-axis// (1,1,1) 



FIG. SI: NV Schematic, (a) Depicts the diamond unit cell 
with corresponding NV axis || (1, 1, 1). The underlying 
lattice of the 2D array can be chosen with xl and yL 
defining the lattice grid (where {xc,yc,Zc) represent the 
Cartesian basis). For simplicity, in the derivations, we 
assume a (1,1, l)-cut diamond, whereby the plane of the cut 
is parallel to the computational plane spanned by a;_L and 
yL- (b) In this particular case, all NV centers and Nitrogen 
impurities sit in the computational plane defined by xl, yL- 
We note that this choice of axes is in slight contrast to the 
main text, where for simplicity x, y is chosen to represent 
the lattice in which all spins lie. Dark spins represent 
Nitrogen impurities while NV centers are represented as 
pairs of green (electronic) and yellow (nuclear) spins (in 
analogy to the main text). Here, in contrast to the sparsely 
occupied plaquettes in the main text, we have chosen to 
occupy all potential sites in the computational plane with 
spins for illustrative purposes. 



sa 2ms time of the Nitrogen impurity is limited by dy- 
namic Jahn- Teller reorientation [S3j [SH [S5] . Thus, nom- 
inally, the total time required for both DSCB mediated 
coherent coupling protocols introduces a significant depo- 
larization error if we consider limitation by the impurity 

time. However, at room temperature, the dynamic 
Jahn- Teller (JT) reorientation is governed by tunneling 
between the four tetrahedrally symmetric axes, and the 
associated rate is given by }S6| 



1/rf = 4* 10^ 



(S4) 



where e = 0.76eV is the experimentally determined ac- 
tivation energy. This exponential dependence on both 
activation energy and temperature suggests that a com- 
bination of a static electric field (- lO^V/cm) [STjlSSllSO] 
and slight cooling (T « 250K) can extend to above 
IDs. Hence, in the discussion of operational errors in 
the main text, we have assumed limitation by T^^ ^ the 
spin-lattice relaxation time of the NV center. 



Interaction between NV Quhit and N Impurity 

In addition to the properties associated with single Ni- 
trogen impurities, we now consider the NV-N interaction 
and show that it is possible to distill an Ising interaction. 



The NV Hamiltonian is given by 

HnV = gefJ-sBzSz — gNfJ-NSzI; 

+ Anv{SzIz + Sxlx 



Syly) 



(S5) 



where A^v — 3.1MHz and Ao = 2.87GHz is the zero- 
field splitting of the singlet [ SIO] . For the purposes of this 
derivation, we will consider the = 0,1 states to be the 
NV electronic qubit and will disregard the nis = — 1 state 
which is far-detuned (Zeeman shifted away due to the 
~ IT magnetic field). The fiip terms between the elec- 
tronic and nuclear spin (5*+/^ + 5^/+) of the NV are 
suppressed due to the zero-field splitting and the large 
magnetic field, since Aq, ge^sBz ^ Any- The dipole- 
dipole interactions between the NV electronic spin and 
an impurity electronic spin will be nearly identical to the 
impurity-impurity dipole coupling which has previously 
been considered. In particular, considering quantization 
along the NV axis, which for simplicity will also be the as- 
sumed Jahn- Teller axis (the result is analogous for other 
JT orientations), and noting that the unit vector pointing 
from the NV to the N-impurity lies in the computational 
plane defined by the normal vector (1,1,1), we find that 
the Hamiltonian is (under the secular approximation in 
analogy to Sec. lA), 



H 



-f^w5f^/f^ + Ao(5f^)2 

^0 ^2, ,2 cNV cN 



NV 



(S6) 



g'^^%srsi 



To keep only the Ising interaction term, we require wnv, 
ujn, and \ujnv — WAr| to all be much greater than the 
strength of the dipole coupling, where uj^v is the coeffi- 
cient of S^^ and ujn is the coefficient of . 



Simulations and Error Optimization 

We consider the optimization of errors in the context of 
the DSCB mediated remote register coupling. First, we 
consider the accumulated infidelity associated with the 
adiabatic SWAP sequence, represented by Eq. 5 in the 
main text, p^J^ « N{padia + pf// + Pdip + Pti + PtD- 
Re-expressing pf^- to capture the form of the associated 
infidelity yields 



ferr 



N 



1 



{n't. 



91 



Ti 



NV 



ts 



rpNV 

^2 



(S7) 

where tss represents the time required to SWAP between 
a single pair of impurities, k' « 8.7kHz, VL^ represents the 
Rabi frequency on the impurities, = lOMHz repre- 
sents the gradient-induced splitting for a magnetic field 



gradient 
(see Sec. 



lO^T/r 



rpNV 



= 250ms, and ^ = 10ms 



lie for details on Ag derivation). Here, we 



describe the error terms in more detail to evince the 
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origin of the infidehty. The first term in pf^^ corre- 
sponds to the non-adiabatic correction resulting from an 
optimized adiabatic ramp profile [Sill IS12[ IS13j . The 
second term represents off-resonant excitations induced 
by The third term corresponds to additional off- 
resonant errors induced by the finite initial splitting as 
well as by dipole-dipole coupling of next-to-nearest neigh- 
bor impurities; since the dipolar interaction strength ~ k' 
while the characteristic energy spacing between impuri- 
ties ^ Vli, the associated error is [ff^ ■ The fourth 
error term corresponds to the depolarization error in- 
duced by the finite NV Ti time, tss/T{^^. The final error 
term corresponds to the infidelity induced by dephasing 

The combined total error can be numerically opti- 
mized. However, before discussing the results of this op- 
timization, we note that the limiting time scale in the sys- 
tem corresponds to ; in particular, error estimates 
and numerical optimization are nearly independent of 
T^^, so long as > 10ms, which is experimentally 

achievable in isotopically purified diamond samples |S15| 
. Furthermore, we note that the current form of pt2 as- 
sumes that only a single echo pulse is applied during an 
individual SWAP. This suggests that the functional form 
of the infidelity can be significantly improved by increas- 
ing the number of echo pulses, which is in effect already 
achieved in our system due to the assumed strong driv- 
ing [ S16[IS17) . Hence, in the remaining error analysis, we 
assume that dephasing-induced errors can be neglected. 
Numerically optimizing the full infidelity in param- 
eters tss and rti (for N = 18), yields fti w 450kHz and 
Ntss ~ 3ms with total error ss 2.6 * 10^^, in good agree- 
ment with the full numerical simulations presented in the 
main text. The numerical simulations account for the ad- 
ditional Jahn- Teller (MW) frequencies and arc obtained 
through numerical integration of the Schrodinger equa- 
tion with the optimized adiabatic SWAP profile. Non- 
unitary errors corresponding to depolarization are then 
added to the unitary errors and the total infidelity is sub- 
sequently optimized. 

Next, we consider the accumulated infidelity associ- 
ated with FFST, represented by Eq. 6 in the main text 

P^//^ ^P^f/^+Pferrm+Pg+PrC^^ (nCglccting 

as discussed above) . Re-expressing p^^^'^ to capture the 
form of the associated infidelity yields 



FFST 



02 



(A,) 



— +( 



k/N ' ^^a' ^ Tf^ ' 



(SB) 



of an NV register. Similarly, the second term also results 
from an off-resonant error and corresponds to undesired 
coupling with off-resonant fermionic modes |S18j . Since 
the coupling strength is characterized by g/VN, while 
the characteristic eigenenergy splitting ~ n/N, such an 

off- resonant error induces an infidelity ( ^^J^ )^. The 
third error term results from the protocol designed to 
achieve controlled coupling g, as elucidated in the Meth- 
ods. Finally, directly analogous to pf^, the final term 
corresponds to the error induced by the operational time, 
tFFST ~ VN/g, which results in depolarization. How- 
ever, it is essential to note that relative to the adiabatic 
sequential SWAP, there is an additional factor of N in 

„FFST 

Pti 

This factor results from the generation of multi-partite 
entanglement via a set of controlled-phase gates during 
a single FFST step; despite such entanglement, remote 
register coupling is readily achieved because a second 
transfer step (corresponding to the return of the quan- 
tum information to the original NV as shown in Fig. 3C 
of the main text) disentangles the quantum information 
from the intermediate dark spin chain data bus |S18j . 
Crucially, the local two-qubit gate which is performed in 
the middle of two FFST steps must commute with the 
controlled-phase multi-partite entanglement. In particu- 
lar, the two-qubit operation depicted in Fig. 3C could be 
a local CP-gate, which would be preserved in the disen- 
tangling step. Combined with single qubit rotations, a 
controlled-phase gate allows for arbitrary two-qubit gates 
and hence universal computation. 

The controlled coupling g is achieved by utilizing the 
NV three-level structure, whereby g ^ k£ (as discussed 
in Methods section of the main text), allowing for the re- 
expression of g in pfj,^^'^ . Similarly, since FFST requires 
tuning to a particular fermionic eigenmode, it is also pos- 
sible to re-express A fa fijv (see Sec. IIB for details on 
fermion tuning). These re-expressions allow for numeri- 
cal optimization in the two parameters VIn and yielding 
Hn ^ 285kHz and ^ 95kHz, with total operational er- 
rors ~ 2.4* 10^2 and total transfer time tpFST = 0.21ms 
(where we have utilized parameters: T-^^ = 250ms and 
N = 7), in good agreement with the full numerical sim- 
ulations presented in the main text. The numerical sim- 
ulations of FFST also account for the additional Jahn- 
Teller frequencies and are obtained through full Hamilto- 
nian diagonalization. Non-unitary errors corresponding 
to depolarization are then added and the total infidelity 
is subsequently optimized. 



where 51 at corresponds to the Rabi frequency applied on 
the impurity, 57 corresponds to the Rabi frequency ap- 
plied on the NV register, k « 12.6kHz, tppsr represents 
the total time required for state transfer, and A corre- 
sponds to the NV detuning as defined in the Methods sec- 
tion in the main text. Here, we describe the error terms in 
more detail to evince the origin of the infidelity. The first 
term in p^^^"^ corresponds to the off-resonant excitation 



Specific Implementation of Architecture 

Here, we offer a specific implementation of the architec- 
tural design principle and discuss the various ingredients 
required to achieve DSCB-mediated coherent coupling 
between spatially separated NV registers. In particular, 
we consider the refocusing of non-nearest neighbor in- 
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FIG. S2: Depicts the design principles governing the 
architecture. Specific impurity implantation within the 
architecture, which is designed to achieve parallel gate operations. 
Plaquettes, outlined in grey, are chosen with y 650nm and 
X ^ 525nm; the individual plaquette-lattice spacing is chosen with 
h 6nm and w ~ 19nm. The distance between all nearest 
neighbor pairs in the horizontal direction is 20nm (-s/h? + w)^ = 
^JXZK)^^+l^v^p2^) and the x spacing between spin D and spin A' 
is given by w/2. We note that the figure is not drawn to scale and 
that all impurity spins sit on the horizontal lattice, while some 
impurity spins sit at midpoints of the vertical lattice. The 
particular lattice spacing is chosen to ensure that the spatial 
separation between individual Nitrogen impurities forming the 
spin chain, as well as between NV registers and Nitrogens does not 
exceed 20nm. The chosen plaquette-lattice spacing implies that 
neighboring plaquettes are coupled by 25 Nitrogen impurities 
(black spins) in the x direction and by Ri 25 Nitrogen impurities 
in the ij direction. The zig-zag pattern of impurity implantation 
in the y direction {w/2 horizontal distance) ensures that all spin 
chain links are of identical length, thereby allowing for parallel 
sequential SWAP operations along y. Additionally, such a pattern 
overcomes the limitation of imperfect NV conversion efficiency by 
allowing each vertical spin chain to span a sufficiently large 
horizontal spacing. Similarly, the saw-tooth pattern of impurity 
implantation in the x direction is designed to allow for echo pulse 
sequences which enable the cancelling of next-to-nearest neighbor 
interactions. In addition to refocusing next-to- nearest neighbor 
interactions, since the impurities forming the horizontal spin 
chain occupy unique rows relative to the vertical spin chain, 
during the adiabatic sequential SWAP, selective spin echoing of 
the horizontal chain impurities ensures that dipole coupling from 
such impurities will not induce operational errors during the 
vertical SWAP sequence. Similarly, selective spin echoing of 
vertical spin chain impurities during FFST can also suppress 
analogous operational errors. Finally, while the architecture is 
designed to showcase the ability to utilize the adiabatic sequential 
SWAP and FFST in differing directions, it is important to note 
that within the specific blueprint above, either implementation of 
the DSCB can be used in both directions. 



teractions along the horizontal spin chain, the control of 
directionality during FFST, and the full super-plaquette 
frequency spectrum. 

Refocusing Next- Nearest- Neighbor Interactions 

The effective Hamiltonian evinced in Eq. 4 of the main 
text has nearest neighbor form; although next-nearest- 
neighbor interactions will represent a correction, we show 
that such interactions can be refocused within the cur- 
rent architecture design, and further, that in principle, 
interactions beyond next-nearest-neighbor can also be re- 
focused. In particular, the horizontal spin chain [N total 
spins) is arranged in a staggered saw-tooth fashion, as 
shown in Fig. S2. Within such an architecture, nearest 
neighbor coupling terms correspond to all pairs of adja- 
cent spins, each separated by ~ 20nm, with correspond- 
ing interaction Hamiltonian 

Hiq — ^ ^ i^n{S\S Q + S^S(j + SqSj^ + S^S^,) -f h.c. 

. 

where the sum runs over all nearest neighbor pairs in 
a given dark spin chain. Thus, next-to-nearest neigh- 
bor terms for each spin correspond to the subsequent 
strongest interaction 

Hnn = ^ knn{S^S^ + S^S];,) ^^^^^ 
+ UNN'iS^S^^, -\- S'^Sg,) + h.c. 

where the prime denotes the next link in the saw-tooth 
chain as shown in Fig. S2. In addition to the impurity 
spins, FFST incorporates the electronic spin of the NV 
register into a mixed spin chain. It is important to note 
that the spin-flip Hamiltonians Hj^ and Hpfj^ are de- 
rived from the secular approximated Ising coupling by 
the application of driving fields as per Eq. 3 in the main 
text. Since each row (1,2,3,4) is separately addressable 
by virtue of the magnetic field gradient (again, applying 
four frequencies for each row to ensure that all JT and 
nuclear spin states are addressed), it is possible to apply a 
spin-echo procedure to refocus the next-nearest-neighbor 
terms. In particular, by flipping the spins in rows 1 and 
2 (Fig. S2) after time Td/2 where Td is a small frac- 
tion of the desired evolution duration, the next-nearest- 
neighbor interactions are refocused since each term con- 
tains spins from only row 1 or 2. However, half of the 
nearest neighbor interactions are also refocused, leav- 
ing effective evolution under the Hamiltonian -ffe//i — 
J2 KN{S^Sg -\- S'^ Sjy) -hh.c. Analogously, by flipping the 
spins in rows 2 and 3, effective evolution under the Hamil- 
tonian Heff2 = f^N{SgS^ + 5*^5^,) +h.c. is achieved, 
again with Hnn refocused. Combining the evolution ac- 
cording to Hgffi and ^^e//2 yields the desired nearest- 
neighbor Hamiltonian with next-to-nearest neighbor in- 
teractions refocused. However, since H^ffi and -ffe//2 
do not commute, it will be necessary to employ piece- 
wise evolution according to the Trotter-Suzuki formal- 
ism |S191 IS201 IS21| . Further refocusing of higher order 
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non-nearest neighbor interactions can also be achieved 
by extending the number of rows corresponding to the 
saw-tooth design; such an extension allows for the isola- 
tion of each specific pair of nearest neighbor interactions, 
thereby achieving the desired nearest-neighbor evolution 
through a Trotter sequence. 

Control of Directionality in Free Fermion State Transfer 

We consider the Hamiltonian presented in Eq. 4 in the 
main text. Under the assumption that g <C '^6 
work perturbatively in eigenstates of the bare Hamilto- 
nian ffo = Y^fji i^{SfS~^-^ + S~Sf_^-^), and consider 
coupling through the perturbation Hamiltonian H' = 
g{Sjjy^S^ + Stjv^SJi + h.c. ). The essence of FFST can 
be understood as the long-range coherent interaction be- 
tween the spin qubits, mediated by a specific collective 
eigenmode of the intermediate spin chain. This mode 
is best understood via Jordan- Wigner (JW) fcrmioniza- 
tion [S2l [S23l [S24] . which allows for the states of an 
XX spin chain to be mapped into the states of a set of 
non-interacting spinless fermions. In this representation, 
the state transfer is achieved by free fermion tunneling. 
By ensuring that the end spin qubits are weakly coupled 
to the intermediate spin chain, it is possible to tune the 
NV registers to achieve resonant tunneling through only 
a single particular fermionic eigenmode jSlSj . Particle- 
hole symmetry of the Hamiltonian implies that the en- 
ergy spectrum of the single fermion manifold is mirror 
symmetric across E = Q. This implies that in the case 
of even N intermediate spin chains (which are uniquely 
present in the proposed architecture), the NV registers 
are always initially off-resonant from all fermionic eigen- 
modes of both left and right spin chains. 

In this context, directionality becomes easily achiev- 
able so long as the left and right spin chains are of differ- 
ing lengths, since the single fermion spectra will then 
be distinct. By ensuring that g is sufficiently small, 
it is possible to tune to and hence be resonant with 
only a single spin chain direction. In particular, the en- 
ergy spectrum of the single fermion manifold is given by 
Ek = 2kcos where k = I,-- - , iV; thus, tuning to 

a particular fermionic eigenmode corresponds to making 
choices of A and VI which ensure that A — Sljv ~ -Efe . The 
specific choice of k corresponds to the particular single 
fermion eigenmode which is being tuned to. Interest- 
ingly, this can allow for control over the speed of quan- 
tum state transfer |S18j . This speed is maximimized for 
k = ^ ±1 in the case of even N chains and for k = 
in odd N chains. Directional control over state transfer 
is achieved by ensuring that only either the left or right 
DSCB is resonant with the NV register. In this scenario, 
the coupling between the NV and the neighboring spin 
chain, which is ofi-resonant is highly suppressed. Assum- 
ing that the two neighboring spin chains are of differing 
lengths A^i and N2, the characteristic energy separation 
between fermionic eigenmodes in the two chains is ap- 



proximately = Thus, by ensuring that 

the register-impurity coupling g is smaller than such an 
energy separation, it is possible to ensure that only single 
directional FFST occurs. Additionally, such an analysis 
suggests that by tuning g, it may be possible to overcome 
coupling-strength disorder induced by imperfect impurity 
implantation |S18| . In particular disorder will cause lo- 
calization, asymmetry of the eigenmodes, and changes in 
the statistics of the eigenenergies. In the case of coupling- 
strength disorder, there exists an extended critical state 
at -E = with a diverging localization length; this ensures 
the existence of an extended eigenmode with a known 
eigenenergy, suggesting that FFST is intrinsically robust 
against couphng-strength disorder |S25[ IS26] . However, 
the existence of an extended mode is not sufficient to en- 
sure state transfer as disorder also enhances off-resonant 
tunneling rates and causes the eigenmode wavefunction 
amplitude to become asymmetric at the two ends of the 
chain. Despite such imperfections, by individually tuning 
the qubit-chain couplings, gie/t and gright, it is possible 
to compensate for eigenmode asymmetry; furthermore, 
sufficiently decreasing the magnitude of the qubit-chain 
coupling ensures that off-resonant tunneling can safely 
be neglected, even in the presence of disorder [S18) . 

Full Frequency Spectrum in Super-plaquettes 

As mentioned in the main text, the frequency spectrum 
of a given super-plaquette requires careful consideration. 
In particular, given a sufficiently large magnetic field gra- 
dient, it will be essential to consider the full range of fre- 
quencies covered throughout a super-plaquette. For ex- 
ample, although the NV zero-field splitting (ZFS) Aq ~ 
3GHz naively protects the register from off-resonant exci- 
tations during impurity manipulation, it is important to 
note that this static ZFS is bridged as rows progress ver- 
tically along the field gradient. Given a gradient of order 
lOOMHz/lOnm (in this section, for simplicity, we assume 
lOnm row separation), a displacement of only 30 rows 
along the gradient direction is sufficient to fully bridge 
the NV ZFS. Furthermore, since each Nitrogen impu- 
rity requires four separate frequencies to be addressed, 
it is essential to ensure that such a congested frequency 
spectrum within a super-plaquette does not render the 
error probability intolerable. As discussed, the four fre- 
quencies corresponding to each Nitrogen impurity result 
from the two possible nuclear spin states and the choice 
of whether the JT axis lies along the NV axis; hence, 
the possible frequencies are ujq ± l/2yl|| ~ uq ± 80MHz 

and Wo ± l/2(4f + ^) ^ ^0 ± 60MHz. To describe 
an adequate framework regarding the frequency spec- 
trum, it is possible to consider a simplified problem, in 
which we ignore the nuclear spin and thus, are left with 
three transitions corresponding to the two NV transitions 
(|0) |±1)) and the Nitrogen impurity transition. For a 
static magnetic field that gives a Nitrogen ESR frequency 
equal to ujQ, the NV transitions correspond to loq + 3GHz 
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and Wo — 3GHz. Thus, we can describe the situation 
according to three base hues at —3, 0, 3GHz, where the 
full frequency spectrum within the super-plaquette re- 
sults from the addition of all possible lines at subsequent 
lOnm intervals. 

One possible solution to avoid frequency overlap is to 
choose the field gradient to be 3 * CMHz/lOnm, where C 
is as yet undetermined. The possible choices of C. corre- 
spond to field gradients which for some n G Z yield that 
after n * lOnm, the field due to the gradient is 3GHz— C. 
Such a choice immediately ensures that (disregarding JT) 
all possible transitions within a super-plaquette are non- 
overlapping. Furthermore, by correctly choosing it is 
possible to ensure that even within the context of ad- 
dressing the four individual Nitrogen impurity transi- 
tions, the separation of any pair of nearest frequencies is 
~ lOMHz. A simple example can be illustrated by con- 
sidering a gradient with strength 150MHz/10nm, where 
the absolute minimum spacing between any two frequen- 
cies in a super-plaquette is Ag = lOMHz. 



Proposed Experimental Realization and Limitations 

NV Implantation and Conversion 

Although we have considered various errors accumulated 
during single- and two-qubit gate operations, an impor- 
tant further consideration is the error associated with the 
imperfect positioning of Nitrogen ion implantations. This 
imperfect implantation leads to errors in both the cou- 
pling strength as well as the individual impurity ESR fre- 
quencies. Crucially, both the vertical adiabatic sequen- 
tial SWAP and the horizontal FFST are robust to fluc- 
tuations in coupling strength. The adiabatic sequential 
SWAP is also robust to errors induced by ESR frequency 
variations. Such variations result in an effective Rabi 



frequency, which alters the start and end point of the op- 
timized ramp profile, without significantly affecting the 
adiabatic passage |S11| . On the other hand, such ESR 
frequency fiuctuations induce an off-resonant error dur- 
ing FFST; however, improvements in the implantation 
precision and utilizing larger D, n can sufficiently suppress 
such errors. 

In addition to implantation errors, it is also important 
to consider the infidelity induced by imperfect NV con- 
version efficiency. In particular, current experiments are 
limited by an optimistic NV conversion efficiency of ap- 
proximately 40%. Within the context of the current ar- 
chitecture design (Fig. S2), each plaquette contains 8 
possible implantation sites that will allow for a stag- 
gered super-plaquette NV register structure, yielding a 
nominal error below 10~^. Thus, each super-plaquette 
will have a filling fraction greater than 99% of functional 
plaquettes and the gate overhead associated with such 
faulty qubits will be negligible [S27j . Since the threshold 
penalty for faulty qubits is expected to be proportional 
to the gate overhead, the error threshold e is expected to 
remain « 1.4% [S27..S28) . 

Furthermore, the zig-zag structure of the vertical impu- 
rity spin chain allows for connectivity between any NV 
lattice site and the lattice site directly above in the corre- 
sponding vertically adjacent plaquette, as well as any of 
the four nearest neighbor lattice sites in either direction 
(in order to account for imperfect conversion efficiency), 
as can be seen in Fig. S2. The implantation region of 
vertically adjacent plaquettes is somewhat limited by the 
particular location of NV registers in the preceding ver- 
tical plaquette; however, such errors can be made negli- 
gible by considering larger plaquette sizes and hence, a 
larger range of connected vertical sites through the im- 
purity spin chain. 
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